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Aanleiding  en  doelstelling  van  het  onderzoek 

Door  het  TNO-FEL  is  een  beeldbewerkingsalgoritme  ontwikkeld  op  basis  van  adaptieve 
kontrastversterking  met  het  oog  op  het  vergroten  van  het  waameembereik  van  camerasystemen 
onder  gevechtsomstandigheden.  Het  is  gebleken  dat  onder  alle  omstandigheden  beeldbewerking 
met  dit  algoritme  tot  een  komfortabeler  waameming  en  vaak  tot  een  vergroting  van  het 
zichtbereik  leidt.  Onder  bepaalde  omstandigheden  kan  het  zicht'  met  een  faktor  1,5  groter 
worden,  bijvoorbeeld  van  4  tot  6  km.  Evaluatie  van  deze  bewerking  als  functie  van  diverse 
atmosferische  omstandigheden  is  van  groot  belang.  Deze  evaluatie  is  onderwerp  van  in  dit  rapport 
beschreven  studie. 

Omschrijving  van  het  onderzoek 

Een  groot  aantal  1 2-bit  beelden  zijn  zo  bewerkt  met  dit  EEL  beeldbewerkingsalgorithme  dat  een 
optimale  presentatie  van  de  door  de  camera  geregistreerde  informatie  werd  verkregen.  De  beelden 
zijn  opgenomen  met  het  eerste  FEL  12-bit  CCD  camerasysteem  vanaf  de  toren  van  het  TNO-FEL. 
De  resultaten  van  de  bewerking  worden  beschreven  en  toegelicht  met  een  aantal  voorbeelden. 
Verbanden  tussen  de  gevonden  beeldbewerkingsparameters  en  de  atmosferische  omstandigheden 
zijn  onderzocht.  Een  duidelijk  veiband  werd  gevonden  tussen  het  'laagste  geregistreerde 
lichtniveau'  in  een  beeld  en  de  Telatieve  vochtigheid  gedeeld  door  het  zicht'.  Het  algoritme  bevat 
een  globale  bewerking,  de  'stretch'  transformatie,  en  een  lokale  adaptieve  kontrastversterking.  Bij 
aanwezigheid  van  slechts  een  beperkte  lichtomvang  (mist,  heiig  weer),  wordt  met  een  'stretch' 
transformatie  (op  de  pixeiwaarden  van  het  geheie  beeld  of  van  een  deelbeeld)  alleen  reeds  een 
aanzienlijke  beeldverbetering  verkregen. 
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Conclusies  en  aanbevelingen 

Bij  toepassing  van  het  algoritme,  vooral  bij  slecht-zichtsituaties,  kunnen  meer  details  via 
bewerking  zichtbaar  worden  gemaakt  en  kan  dus  een  vergroting  van  het  zichtbereik  gerealiseerd 
worden.  De  mate  van  zichtverbetering  hangt,  behalve  van  de  camera  performance,  in  de  praktijk 
ook  af  van  de  aard  van  de  scene,  van  aard  en  mate  van  het  kontrastverlies  in  de  atmosfeer  en 
tenslotte  ook  van  de  waamemer.  Voor  slecht-zichtsituaties  (mist  en  heiig  weer)  kan  een 
verbetering  van  het  'zicht'  met  rond  een  faktor  twee  verwacht  worden.  Een  betrouwbare 
kwantitatieve  bepaling  vereist  zeer  geconditioneerde  experimenten. 

De  met  deze  studie  verkregen  inzichten  zijn  benut  voor  andere  CCD  projekten.  Een  12-bits 
camerasysteem  met  PC  besturing  en  met  een  Peltier  gekoelde  camerakop  werd  inmiddels 
gebouwd  om  opnamen  op  willekeurige  lokaties  te  kunnen  maken.  Hiermee  werd  bijvoorbeeld  een 
groot  aantal  opnamen  gemaakt  tijdens  de  'BEST-TWO'  battle  field  trial. 

Voor  toepassingen  is  een  real  time  working  van  de  beeldverbetering  van  groot  belang.  Een  real 
time  hardware-uitvoering  van  een  variant  op  het  algoritme  werd  inmiddels  gerealiseerd  in  de 
vorm  van  een  real  time  videoprocessor.  In  het  kader  van  een  nationaal  technolc'gieprojekt  werd 
eveneens  een  real  time  hardware-uitvoering  van  het  algoritme  gerealiseerd  tezamen  met  een 
aantal  real  time  werkende  filtermodules. 

De  vergroting  van  herkenningskansen  kan  op  een  andere  manier  afhangen  van  de  bewerking  dan 
de  vergroting  van  het  zichtbereik.  Onderzoek  naar  de  vergroting  van  herkenningskansen  van 
diverse  militaire  objekten,  door  toepassing  van  beeldbewerking,  wordt  aanbevolen.  Het  is  van 
belang  het  nut  van  beeldbewerking  verder  te  kwantificeren  en  te  evalueren  onder  diverse 
omstandigheden. 
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ABSTRACT 

A  survey  is  given  of  processing  results  of  12-bit  pictures  with  the  PEL  image  processing 
algorithm.  The  results  are  illustrated  with  examples.  The  pictures  have  been  taken  with  the  first 
PEL  12-bit  CCD  camera  system  under  various  atmospheric  conditions.  Some  relations  between 
processing  parameters  and  atmospheric  conditions  are  analysed.  An  indication  of  a  quantitative 
visibility  range  improvement  is  given. 

SAMENVATTING 

De  resultaten  \an  de  bewerking  van  12-bit  beelden  met  het  PEL  beeldbewerkingsalgoritme 
worden  beschreven  en  toegelicht  met  een  aantal  voorbeelden.  De  beelden  zijn  opgenomen  onder 
diverse  weersomstandigheden  met  het  eerste  PEL  12-bit  CCD  camerasysteem  vanaf  de  toren  van 
het  TNO-PEL. 

Enkele  relaties  tussen  beeldbewerkingsparameters  en  atmosferische  omstandigheden  zijn 
onderzocht.  Voorde  zichtverbetering  als  resultaat  van  de  beeldbewerking  wordt  een  kwantitatieve 
indicatie  gegeven. 
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1  INTRODUCTION 

The  study  to  be  described  in  this  report  is  carried  out  to  evaluate  a  sophisticated  TV-perception 
system  based  on  the  first  prototype  of  a  12-bit  CCD-camera  [1].  The  main  objectives  are:  to 
improve  the  visibility  range  and  apparently  reduce  the  loss  of  contrast  in  the  atmosphere.  At  the 
Physics  and  Electronics  Laboratory  TNO  a  CCD-camera  system  has  been  developed  with  a  12-bit 
Intrascene  Dynamic  Range.  Special  attention  is  paid  to  the  analog  electronic  read-out  circuits  and 
to  the  cooling  of  the  sensor  to  obtain  the  necessary  signal-to-noise  ratio  1 1  ]. 

To  present  all  the  information  of  a  12-bit  image  on  a  monitor,  image  processing  is  necessary.  A 
special  algorithm  has  been  developed,  that  transforms  the  dynamic  range  of  the  images  to  the 
optimum  range  of  the  monitor.  Simultaneously  the  contrasts  are  enhanced  such,  that  all  the 
recorded  information  becomes  within  the  perception  range.  A  reduction  of  the  dynamic  range  also 
results  in  a  reduction  of  the  contrasts.  The  digital  image  processing  therefore  must  do  anyhow  a 
contrast  enhancement.  It  is  also  possible  to  restore  the  loss  of  contrast  in  the  atmosphere  with  this 
contrast  enhancement. 

This  report  mainly  de.scribes  the  evaluation  of  the  image  processing  algorithm  as  function  of 
various  atmospheric  conditions.  In  the  considered  period  many  pictures  have  been  taken  w  ith  the 
TNO-camera  system  and  processed  with  the  algorithm.  A  data-base  of  images  has  been  set  up  and 
this  data-base  has  been  analysed  with  stati.stical  methods  to  retrieve  possible  relationships  with 
meteorological  parameters.  Also  a  general  parameter  setting  of  the  algorithm  is  determined,  such 
that  this  set  can  be  applied  on  images  taken  under  different  conditions.  This  was  done  in  view  of 
possible  future  automation.  A  first  indication  of  a  quantitative  visibility-improvement  is  given.  In 
chapter  two  the  digital  image  processing  algorithm  will  be  discussed  briefly.  For  a  more  extended 
paper  on  this  subject  is  referred  to  the  report  FEL-1988-63  [2]  and  de  Vries  [3].  In  chapter  three  a 
description  of  the  image  data-base  will  be  given  and  the  image  acquisition  will  be  discussed.  A 
brief  introduction  to  the  theory  of  visibility  in  the  atmosphere  will  be  given  in  chapter  four.  It  will 
be  shown  that  carefully  defined  experiments  are  necessary  to  quantify  the  'visibility  range 
improvement'.  The  results  will  be  discussed  in  chapter  five  and  finally  in  chapter  six  some 
conclusions  will  be  given.  This  study  has  been  carried  out  mainly  under  the  assignment  number 
A87/KL/I59.  A  critical  analysis  of  the  re.sults  and  a  final  update  of  the  report  has  been  carried  out 
under  the  assignment  number  A9()/KL/675. 
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DIGITAL  IMAGE  PROCESSING 


2.1  Introduction 

In  this  chapter  the  algorithm  used  to  process  the  images  will  be  described  shortly.  For  a  more 
extended  paper  on  this  subject  is  referred  to  the  FEL-report  FEL-1988-63  [2].  This  algorithm  for 
digital  image  processing  has  been  developed  at  the  Physics  and  Electronics  Laboratory  TNO. 
Also  a  few  pre  process  and  post  process  procedures  will  be  described.  Finally,  an  optimum 
procedure  for  processing  the  images  with  the  same  set  of  processing  parameters  will  be  given. 

2.2  The  image  processing  algorithm 

The  purpose  of  the  iiiiage  processing  is  to  bring  the  information  recorded  with  the  12-bit  CCD- 
camera  into  the  perception  domain  of  a  monitor.  It  is  necessarj'  therefore,  to  reduce  the  dynamic 
range  from  12  bits  to  8  bits;  it  is  the  common  digital  input  for  standard  TV  applications.  This 
transformation  w  ill  result  in  a  decrease  of  contrast.  The  task  of  the  image  processing  therefore,  is 
to  enhance  the  small  contrasts,  which  can  be  due  to  the  loss  of  contrast  in  the  atmosphere  as  well 
due  to  the  dynamic  range  reduction. 

Of  course  commercially  available  image  processing  systems  offer  a  score  of  relative  simple 
techniques  to  enhance  a  digital  image.  A  common  disadvantage  of  these  techniques  is  that  often 
only  a  part  of  the  image  is  enhanced  at  the  cost  of  another  part.  The  here  used  technique  operates 
on  the  full  image  and  carries  out  the  enhancement  independent  of  the  local  brightness. 

The  range  transformation  is  calculated  according  the  equation  <2.1>: 

Y  =  <2.1> 

X  is  the  grey  level  in  the  input  domain  and  Y  in  the  output  domain. 

The  technique  has  the  potential  of  real-time  u.se  because  each  pixel  is  processed  separately  and 
the  processing  does  not  require  any  kind  of  transform  and  is  not  recursive.  The  dynamic  range  is 
transformed  in  a  non-linear  way,  such  that  the  resulting  contrast  reduction  is  independent  of  the 
local  grey-level. 

The  equivalent  implementation  of  this  tran.sformation  is  .shown  in  figure  2.1 .  This  implementation 
has  the  attractive  property  that  the  local  difference  in  the  logarithmic  domain  is  equal  to  the 
contrast  in  the  input  image,  so  the  contrasts  can  directly  be  manipulated  then. 
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■  igure  2. 1 :  The  implementation  of  the  dynamic  range  transformation. 

The  detail  contrast  enhancement  is  carried  out  by  means  of  a  high  pass-filtering  technique  that 
makes  use  of  only  local  statistics. 

A  local  difference  image  is  obtained  in  the  logarithmic  domain  that  can  be  considered  as  a 
"contrast  image".  The  local  difference  image  is  processed  in  an  adaptive  way  (the  small  contrasts 
are  enhanced  more  than  the  large  ones  to  prevent  artefacts). 

An  implementation  of  the  filter  is  shown  in  figure  2.2. 


z 

+ 

adaptive  local 
difference 

adaptive 

lowpass  filler 

1 

adaptive  local 

mean 


Figure  2.2:  The  implementation  of  the  filtering  technique. 


A  more  extended  description  of  this  filter  can  be  found  in  FEL-report  1988-63  [2].  A  combination 
of  the  figures  2.1  and  2.2  gives  the  complete  processor,  which  is  shown  in  figure  2.3. 
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Figure  2.3:  A  block  diagram  of  the  digital  image  processor  In  this  diagram  U'  is  the  filter  dimension.  \  „ 

is  a  parameter  to  adjust  the  adaptation,  y^.  is  the  contrast  parameter  and  yj  and  k^,^„  are  the 
dynamic  range  reduction  parameters. 

The  optimum  value  of  the  adaptation  parameter  i.s  determined  from  the  local  v  ariance  range  that  is 
present  in  the  image.  .A  thre.shold  value  can  be  defined  from  that  range  beyond  which  'he  level 
of  contrast  enhancement  depends  on  the  local  variance. 

Another  advantage  of  this  algorithm  is  its  ability  to  enhance  the  high  frequencies  and  to  reduce 
the  low  frequencies.  It  means  that  detail  contrast  is  enhanced,  but  slow  ly  vary  ing  differences  in 
luminance  will  be  eliminated.  These  variations  in  luminance  can  result  from  the  atmosphere  or 
can  be  due  to  spatially  varying  scattering  or  transmission  in  the  lens  and/or  the  shutter  of  the 
camera. 

2.3  Pre  processing  and  Post  processing 

Before  images  with  the  above  described  algorithm  are  processed,  a  few  preprocessing  methods 
can  be  applied.  Some  very  common  methods  will  be  described  here,  such  as  dark  field 
subtraction,  flat  field  correction,  median  filtering  and  histogram  modification.  Post  processing  eg 
can  be  the  gamma  correction  of  a  TV-monitor  and  correction  of  a  possible  geometrical  distortion; 
this  latter  eg  can  arise  from  differences  in  the  dimensions  of  the  pixels  and  those  of  the  image 
section  of  the  CCD-sensor. 

2.3. 1  Dark  field  subtraction 

Because  the  effect  of  the  temperature  is  not  uniform  across  the  surface  of  the  CCD,  the  thermal 
generation  of  electrons  also  is  a  function  of  the  position.  This  may  cause  a  noisy  picture.  By 
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means  of  dark  field  subtraction,  which  is  a  ver>'  common  operation,  the  rms  value  of  the  thermal 
noise  and  fixed  pattern  noise  (due  to  fixed  differences  in  the  pixels )  can  be  removed. 

The  dark  field  image  can  be  recorded  with  a  closed  shutter  and  with  maintaining  all  the  other 
parameler  settings  of  the  camera. 


Y  =  0 

Y  =  X-Xj3r^.  otherwise 


<2.2> 


.Although  the  mean  level  of  the  thermal  noise  w-ill  be  reduced,  the  standard  deviation  w  ill  increase 
by  a  factor  v2. 

Some  images  have  been  preprocessed  by  this  method  but  no  clear  improvements  could  be 
observed.  Obviously,  the  mean  level  of  the  thermal  noise  was  already  low  enough.  This  must  be 
due  to  the  cooling  of  the  sensor  (during  the  recording  sessions)  and  to  the  sophisticated  read  out 
and  signal  processing  electronics. 


2..V2  Flat  field  correction 

The  photo  sites  of  the  CCD-sensor  vary  in  sensitivity  in  a  manner  similar  to  random  noise. 
However,  the  individual  photo  sites  each  have  fixed  sensitivity,  so  the.se  differences  in  sensitiv  ity 
can  easily  be  corrected,  eg  bv  recording  a  picture  with  a  uniform  exposure  of  the  CCD-sensor, 
After  dark  field  subtraction  such  an  image  contains  variations,  which  correspond  with  the  actual 
variations  in  sensitivity  of  the  photo  sites  and  with  the  variations  in  exposure  (caused  eg  by  spatial 
variations  in  transmission  of  lens  and/or  .shutter).  With  the  operation  given  in  formula  2..^.  a  fiat 
field  correction  can  be  carried  out  then. 

Y  =  0  if  Xfiat  =  0 

<2.3> 

Y  =  X*(Xavg/xnat)  if  Xfiat  >0 

Xfiat  is  the  individual  pixel  value  in  the  flat-field  image  and  Xavg  is  the  averaged  pixel  value  in 
the  full  fiat-field  image. 

Improvements  will  be  observed  only  if  the  variations  contain  high  frequency  components, 
because  the  low  .spatial  frequencies  will  be  suppressed  already  by  the  image-processing 
algorithm.  In  processing  the  hereafter  mentioned  image-database,  no  significant  improvements 
have  been  observed  with  this  flat  field  correction. 
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2.3.3  Median  filtering 

This  technique  can  be  used  to  remove  the  effects  of  defective  pixels  within  the  CCD-sensor.  The 
values  of  the  defective  pixels  are  exchanged  by  the  median  of  the  surrounding  pixels.  However, 
this  can  be  done  only  if  the  positions  of  the  defective  pixels  are  known. 

2.3.4  Histogram  modification 

The  optimum  illumination  of  the  sensor  with  respect  to  its  S/N  ratio  is  an  exposure  of  the  sensor 
up  to  its  saturation  level  (see  chapter  four).  All  the  pixel  values  can  be  multiplied  with  a  factor  up 
to  the  maximum  value  of  4095.  in  the  case  when  the  sensor  is  underexposed.  Processing  might 
give  better  results  then.  The  effect  of  exposure  is  multiplicative.  In  the  future,  the  histogram  of  the 
image-data  can  possibly  be  used  for  an  optimum  real  time  exposure  of  the  sensor.  In  the  case 
when  the  luminance  range  of  the  image  is  very  limited,  then  better  processing  results  can  be 
obtained  by  carving  out  first  a  range  transformation.  Sometimes  .at  dense  fog.  the  luminance 
range  is  not  more  than  I  or  2  bit  as  a  result  of  scattering  of  light  in  the  atmosphere. 

2.3.5  Correction  geometrical  distortion 

The  image  field  is  determined  only  by  the  focal  length  of  the  lens  and  the  dimensions  of  the 
sensor  (image  section).  Opaque  parts  within  this  image  section,  as  there  are  the  vertical  transport 
CCU-registers  in  an  interline  transfer  CCD-sensor  (see  fig.2.4),  do  not  change  the  field  of  view. 
Because  of  these  opaque  parts,  the  number  of  pixels  in  both  directions  is  not  proportional  to  the 
corresponding  dimensions  of  the  sensor  and  a  geometrical  distortion  will  appear  when  the  digital 
data  of  the  image  are  displayed  directly  on  a  monitor.  The  distance  between  two  adjacent  pixels 
on  the  monitor  is  the  same  in  horizontal  and  vertical  direction  then. 

For  the  here  used  Fairchild  sensor  222  the  horizontal  and  vertical  dimensions  are  respectively 
1 1 .4  and  8.8  mm  and  the  corresponding  number  of  pixels  respectively  are  380  and  488.  Because 
of  these  ratios,  the  image  field  dimension  in  the  horizontal  direction  must  be  multiplied  by  a 
factor  1 .66. 

We  have  applied  only  a  correction  factor  of  about  1.3,  to  be  able  to  display  4  pictures  on  a 
1024x1024  screen.  The  photographs  given  in  chapter  5  have  these  dimensions  therefore. 
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Figure  2.4:  Fairchild  222  pixel  dimensions  (horizontal  x  vertical  =  30  x  1 8  |im). 
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3  THE  IMAGE  DATA  BASE 

3.1  Introduction 

The  here  described  database  of  images  has  been  set  up  with  the  (first)  experimental  CCD-camera 
system  of  the  Physics  and  Electronics  Laboratory  TNO  [1].  An  artist's  impression  of  this  system 
is  given  in  fig.3.1.  All  the  pictures  have  been  taken  under  the  various  atmospheric  conditions 
occurring  in  1988  and  1989.  The  pictures  were  all  taken  from  the  1 1'*’  floor  of  the  EEL  tower. 
This  database  is  mainly  used  for  the  evaluation  of  the  described  image  processing  algorithm.  The 
image  processing  algorithm  is  extensively  described  in  reference  |2]  and  shortly  reviewed  in  the 
previous  chapter.  First  the  image  acquisition  will  be  discus.sed  and  then  a  survey  of  the  available 
images  will  be  given. 

3.2  Image  acquisition 

All  the  pictures  have  been  recorded  with  a  CCD-camera  system  based  on  the  Fairchild  CCD  222 
sensor.  The  CCD-sensor  is  a  rectangular  array  of  384*488  photosites  (see  figure  2.4).  The 
performance  of  a  12-bit  intrascene  dynamic  range  ( 1  ]  could  be  realized  by  cooling  the  chip  down 
to  -3()°C;  by  compensating  the  dc-level  drift  of  the  video  output  and  by  applying  special  designed 
analog  electronics,  like  'correlated  double  sampling'. 

All  the  pictures  ha\e  been  taken  with  a  75  mm  focal  length  and  a  fixed  exposure  time  of 
30  msecond.s.  This  was  the  shortest  shutter  time  to  get  a  nearly  homogeneous  exposure  of  the 
sensor.  Part  of  the  images  was  recorded  with  a  photopic  filter  (or  a  2  mm  BG18  glass)  in  front  of 
the  lens  to  simulate  the  eye  sensitivity.  This  filter  can  advantageously  be  used  when  the  visibility 
is  good,  especially  in  sunlit  scenes.  Without  this  filter  the  near  infrared  radiation,  which  is 
strongly  reflected  by  vegetation,  can  reduce  the  local  contrasts  because  of  local  saturation.  Under 
conditions  of  limited  visibility,  due  to  atmospheric  haze,  the  filter  should  not  be  used  because  of 
the  lower  extinction  coefficient  for  the  longer  wavelengths  of  near  infrared  radiation  (see  chapter 
3).  For  getting  an  improved  S/N  ratio,  up  to  16  images  could  be  added  with  the  modular  camera 
system  1 1 J. 


TNO  report 


TNO  report 


Page 

16 


To  get  an  optimum  recording  of  images,  it  is  important  to  expose  the  sensor  such  that  those  sensor 
elements,  corresponding  with  the  brightest  parts  in  the  scene,  are  illuminated  just  till  its  saturation 
level.  This  is  because  of  the  Poisson  distribution  nature  of  the  photon  noise,  which  implies  a  S/N 
ratio  proportional  to  Vn,  n  being  the  number  of  generated  photoelectrons.  It  results  in  the  best 
local  signal  to  noise  ratio  without  contrast  loss  in  the  brighter  parts.  Furthermore,  the  maximum 
luminance  range  is  recorded  then,  which  is  important  especially  in  scenes  with  a  large  dynamic 
range. 

3.3  The  image  database 

All  the  recorded  images  are  pictures  taken  from  the  surroundings  of  the  laboratory.  The 
atmospheric  conditions,  under  which  the  images  have  been  recorded,  also  have  been  recorded. 
With  this  information,  we  have  been  able  to  relate  a  number  of  atmospheric  conditions  with  the 
parameters  of  the  FEL  image  processing  algorithm.  The  following  meteorological  quantities  have 
been  noted: 

The  luminance,  the  temperature,  the  relative  humidity,  the  wind  speed  and  wind  direction,  the 
cloudiness,  the  type  of  weather  and  of  course  the  visibility. 

The  visibility  is  not  always  unambiguously  defined,  as  will  be  discussed  in  chapter  4.  The  'eye 
visibility'  of  the  person  (who  was  always  doing  the  recordings)  has  been  noted  together  with  the 
visibility  measured  with  an  AEG  point  visibility  meter.  The  location  of  this  PVM,  however,  was 
at  the  seaside  boundary  of  the  town  near  the  laboratory,  while  most  of  the  pictures  have  been 
taken  in  the  opposite  direction  (the  so  called'  city  views).  Furthermore,  the  range  of  the  visibility 
that  can  be  reliably  measured  with  such  an  instrument,  is  limited  from  0.5  to  about  1 5  km. 
Recordings  have  been  made  for  visibilities  from  0.1  km  up  to  30  km.  In  the  scope  of  a  possible 
application  in  a  military  environment,  atmospheric  conditions  with  a  visibility  less  than  about  5 
km  are  the  most  interesting.  Nevertheless,  evaluation  of  the  system  for  visibilities  up  to  30  km 
might  also  be  important  because  'facing  light’  often  occurs  with  good  visibility.  Figure  3.2  shows 
a  histogram  distribution  of  the  recorded  images  with  the  corresponding  range  of  the  visibility. 
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Image  distribution  versus  visibility 


Visibility  ( km ) 

Figure  3.2:  Image  distribution  for  the  various  ranges  of  the  visibility. 

In  most  cases  a  considerable  offset  in  the  greylevel  range  was  present.  This  might  partly  be  due  to 
the  location  from  which  the  pictures  have  been  taken.  No  nearby  objects  were  present,  so  in  most 
cases  this  results  in  a  contribution  of  atmospheric  scattering  to  all  parts  of  the  scene.  It  might  also 
partly  be  due  to  scattering  at  optical  parts  in  the  camera  head,  especially  when  scenes  are  recorded 
with  a  large  luminance  range  [9,1 1  ]. 

A  PC-controlled  portable  version  of  the  CCD-camera  system  has  been  built  since  then.  With  this 
system,  pictures  can  be  taken  from  any  location,  such  that  more  different  scenes  can  be  (and  has 
been)  evaluated.  In  appendix  A,  a  survey  is  given  of  all  the  interesting  images,  which  have  been 
recorded  in  the  here  considered  period  (1988-1989).  All  these  pictures  have  been  processed  with 
the  FEL-algorithm  to  evaluate  1 2-bit  pictures  of  natural  outdoor  scenes. 
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4  VISIBILITY  IN  THE  ATMOSPHERE 

The  purpose  of  this  chapter  is  to  give  a  short  description  of  some  problems  connected  with 
atmospheric  properties.  One  of  the  goals  of  this  project  was  to  realize  a  camera  system  with  a 
visibility  range  larger  than  that  of  commercially  available  systems.  An  estimation  of  the  visibility 
range  improvement  will  be  given  and  an  experiment  to  verify  the  performance  of  the  camera  will 
be  described. 

4. 1  Theory 

Several  aspects  of  the  visibility  will  be  discussed  briefly.  The  most  important  are  the  optical 
properties  of  the  atmosphere,  the  amount  and  distribution  of  light,  the  characteristics  of  the 
objects  within  the  field  of  view  and  the  properties  of  the  human  eye.  Problems  connected  with 
visibility'  are  rather  complex.  Most  of  the  information,  which  is  acquired  through  our  sense  of 
vision,  depends  on  the  perception  of  differences  in  luminance  and  to  a  lesser  extent  on 
chromatics.  The  contrast  (in  luminance)  in  case  of  an  isolated  object  surrounded  by  a  uniform  and 
rather  extended  background  is  given  by  the  equation: 

C  =  (L  -  L  )  /  L'  <4.l> 

L  is  the  luminance  of  the  object  and  L’  is  the  luminance  of  the  background.  The  range  of  C  runs 
from  -1  for  an  ideal  black  object  to  o®  (eg  a  bright  source  against  a  dark  background).  Due  to  the 
extinction  of  light  in  the  atmosphere  [5],  the  contrast  decreases  with  increasing  distance,  resulting 
in  a  decrease  of  the  visibility.  The  theory  of  Koschmieder  [5]  describes  the  alteration  of  contrast 
due  to  the  path  through  the  atmosphere.  For  the  visual  range  of  an  object  in  the  horizontal 
direction  (seen  against  the  horizon  sky),  the  apparent  luminance  of  the  object  is  given  by  equation 
4.2; 

L  =  L(,  e'®*^  +  Lh  ( 1  -e'®*^)  <4.2> 

Here  is  the  luminance  of  the  object  observed  nearby,  o  is  the  extinction  coefficient  and  L^,  is 
the  luminance  of  the  horizon.  It  is  assumed  that  the  atmosphere  in  the  horizontal  direction  is 
homogeneous.  On  the  same  conditions  the  contrast  reduction  of  the  object  is  given  by: 


C,  =  Qe-OR 


<4.3> 
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Now  Cp  is  the  initial  contrast  of  the  object  (nearby  view).  The  studies  of  Mie  [5]  provide  a  good 
theoretical  basis  for  the  interaction  of  light  and  particles  in  the  atmosphere,  which  determines  the 
extinction  coefficient.  In  case  of  monodisperse  aerosols  with  no  absorption,  relation  4.4  is  valid; 

a  =  N  7C  a-  k  <4.4> 

N  is  the  number  of  particles  per  cm^,  a  is  the  radius  of  the  particles  and  k  is  the  so  called 
scattering  area  ratio.  When  the  relative  humidity  for  instance  is  high,  the  radius  of  the  particles 
will  be  large  also,  in  most  cases,  and  the  contrast  reduction  will  be  large  then  due  to  a  large 
extinction  and  scattering.  The  extinction  coefficient  also  depends  on  the  wavelength.  Various 
formulas  have  been  proposed  and  they  share  the  property  that  o  increases  with  decreasing 
wavelength.  Rain  should  be  treated  as  a  special  case,  because  the  extinction  is  proportional  to  the 
number  of  drops,  falling  per  unit  area  per  second.  Therefore,  this  parameter  should  be  recorded 
during  the  measuring  sessions.  With  increasing  distance  the  contrast  will  decrease  until  a  certain 
threshold  value  is  reached  at  which  the  object  no  longer  can  be  seen.  Definition  of  this  threshold 
value  might  be  rather  complex.  It  depends  on  many  parameters  and  it  does  not  have  an 
unambiguous  value.  It  depends  on  the  properties  of  the  eye,  but  also  on  a  few  nonmeasurable 
influences,  like  the  physical  health,  the  mental  attitude  and  the  alertne.ss  of  the  ob.server.  There  are 
three  types  of  relevant  information  in  the  scene:  the  contrasts,  the  background  luminance  and  the 
dimensions  of  the  objects,  at  which  the  observer  is  looking  at.  Figure  4.1  gives  the  thresholds 
related  to  these  parameters.  They  resulted  from  experiments,  which  were  organized  during  the 
Second  World  War  at  the  Tiffany  Foundation  in  the  United  States  [5]. 


visual  angle  threshold  (mRads) 
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Figure  4.1:  Thresholds  of  brightness-contrast  and  visual  angle  (mrads)  for  50  percent  detection  (after 
Blackwell)  at  each  of  the  given  background  luminances.  The  stimuli  were  brighter  than  the 
background  |5|. 

The  threshold  dimension  of  the  objects  is  given  by  the  visual  angle  in  the  vertical  direction.  The 
contrast  thresholds  are  given  only  for  sources  brighter  than  the  background. 

It  will  be  clear  that  it  is  difficult  to  define  the  visual  range  for  a  given  situation.  For  most  purposes 
it  is  accepted  to  u.se  the  so  called  meteorological  visibility,  which  is  defined  as  the  distance  along 
which  the  transmission  due  to  the  atmosphere  is  decreased  to  291-.  It  is  the  distance  at  which  an 
initial  contrast  Co  =  I(X)9f  is  reduced  to  a  contrast  Cr  =  2%. 

It  has  been  mentioned  before,  that  the  present  investigation  is  not  focussed  on  the  visibility  range 
improvement  related  to  the  eye  visibility,  but  to  the  visibility  range  of  a  conventional  eight  bits 
TV-perception  system.  In  view  of  the  visibility  range  improvement,  it  is  useful  to  know  the 
smallest  contrast  that  can  be  recorded  by  the  CCD-camera.  Although  monitor  inputs  accept  eight 
bit,  only  about  30  to  60  different  grey  levels  can  be  distinguished,  due  to  the  limited  dynamic 
range  of  the  monitor  screen  and  the  human  observer.  It  means  that  a  contrast  becomes  perceptible, 
only  when  a  proper  contrast  enhancement  is  carried  out  by  the  image  processing. 


I 


TNO  report 


Page 

21 


When  a  (threshold)  contrast  Cr  before  processing,  given  by  formula  4.3.,  with  0  =  0]  and  R  =  R|, 
will  increase  to  Cr,  due  to  the  processing,  then  this  increase  can  be  translated  in  an  apparent 
decrease  of  a  to  02.  A  new  apparent  visibility  range  R2  (with  a  same  final  contrast  Cr)  can  be 
calculated  then  from  this  a2.  Now  formula  4.5  is  valid: 

0]Ri  =  O2R2  <4.5> 

With  formulas  4.3  and  4.5  a  'visibility  range'  improvement  R2  /  R]  can  be  derived  with  the 
variables  Cr,  which  is  the  contrast  before  processing,  Cr',  which  is  the  contrast  after  processing 
and  Co.  being  the  initial  contrast. 

ln(Cr/Co) 

Ri/R,  = -  <4.6> 

ln(Cr'/Cr)  +  ln(Cr/Co) 

In  figure  4.2  a  graph  of  this  'visibility  range'  improvement  is  given  as  a  function  of  the  contrast 
improvement'  Cr'/Cr  for  three  values  of  the  'contrast  reduction'  factor  Cr/Co.  due  to  the 
atmosphere. 


visibility-range  improvement 


Cr/Co=.02  --*•  Cr/Co=.05 - Cr/Co=.01 


Figure  4.2: 


The  visibility  range'  improvement  R2/R1  as  a  function  of  the  contrast  enhancement'  Cr'/Cr 
for  three  values  of  the  atmospheric  contrast  reduction'  factor  Cr/Co. 
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The  maximum  contrast  enhancement,  which  can  be  applied,  depends  on  the  signal  to  noise  ratio 
of  the  sensor  and  on  the  recording  conditions.  The  differences  between  the  three  curves  indicate 
the  dependance  of  the  'visibility  range  improvement'  on  the  perception  of  the  minimum  contrast 
before  processing. 

To  arrive  at  a  reliable  validation  of  the  'visibility  range  improvement'  under  adverse  atmospheric 
conditions,  an  experiment  will  be  necessary,  which  must  be  carried  c  under  well  defined 
conditions.  Recordings  should  be  made  from  black  and  white  screens  witt,  ..  known  reflection 
coefficient.  These  screens  must  be  placed  in  a  straight  line  at  known  distances  and  their 
dimensions  have  to  increase  with  increasing  distance,  such  that  the  viewing  angle  is  the  same  for 
all  the  screens.  Variable  contrasts  of  the  .screens  can  be  realized,  eg  by  projecting  black  and  white 
lines  and  lowering  their  contrasts  by  an  extra  screen  illumination  with  a  calibrated  source.  When 
the  screens  are  selected,  which  have  the  same  contrasts  before  and  after  processing,  then  a  direct 
visibility  range  improvement'  results  for  the  corresponding  operational  conditions, 

Recordings  could  be  made  with  our  12-bit  system  and  also  with  an  8-bit  system,  to  evaluate  the 
image  processing  as  function  of  the  camera  performance. 
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3  RESULTS 

In  this  chapter  the  results  of  the  image  processing  and  the  statistical  analysis  of  the  image 
database  will  be  discussed. 

5. 1  Evaluation  of  the  image  processing 

It  has  already  been  mentioned  before,  that  dark  image  subtraction  did  not  have  a  significant  effect. 
Apparently,  the  applied  cooling  of  the  chip  down  to  -30°C  has  sufficiently  lowered  the  mean  level 
of  thermal  noise.  Also  the  fixed  pattern  noise  shows  to  be  already  low  enough  and.  due  to  the 
sophisticated  electronics  [  I  ],  the  system  noise  level  as  well.  A  flat  field  correction  was  carried  out 
fora  few  images,  but  al.so  did  not  have  much  effect.  Obviously,  the  sensitivity  of  the  photosites  of 
the  used  sensor  is  sufficiently  homogeneous.  However,  in  images  recorded  under  conditions  of 
heavy  fog.  a  characteristic  vertical  oriented  pattern  becomes  perceivable  after  processing  (figure 
.3.4).  This  effect  may  be  caused  by  inhomogeneous  illumination,  possibly  in  combination  with  a 
digital  round  off  Due  to  the  limited  speed  of  the  shutter  blades,  the  brightness  in  a  recorded 
image  is  highest  in  the  middle  and  is  decreasing  towards  the  edges.  Every  image  line  is  formed  by 
an  array  of  discrete  photosites.  In  the  horizontal  direction,  there  is  an  opaque  part  of  12  pm 
between  the  succeeding  .sensitive  parts  of  the  photosites  (see  figure  2.4).  Then  the  variation  of  the 
brightness  along  a  horizontal  image  line  is  not  a  continuous  function.  After  processing, 
rectangular  areas  with  different  brightness  appear.  The  main  spatial  frequency  in  this  structure  is 
lower  than  that  corresponding  with  the  pixel  pitch.  Pre  processing,  based  on  histogram 
modification  techniques,  gave  better  results  in  some  cases.  When  the  dc-level  of  the  brightness 
distribution  in  the  image  was  valued  more  than  about  20(X)  (half  the  saturation  level),  a  better 
result  was  obtained  by  first  applying  a  linear  transformation  of  the  recorded  grey  level  range  to 
the  full  range.  Such  a  high  dc-level  can  arise  due  to  strong  backscattering  at  particles  in  the 
atmosphere.  In  general,  the  range  transformation  as  given  by  formula  2.1  is  not  optimal  when 
only  quite  small  grey  level  ranges  have  been  recorded.  When  a  linear  or  a  logarithmic  range 
transformation  is  applied,  before  the  FEL-image  processing  algorithm  is  carried  out,  then  a  better 
distribution  of  the  greylevels  is  obtained. 

For  a  number  of  images  the  maximum  grey  level  value  was  less  than  4095.  Then  it  showed  to  be 
advantageous  for  optimum  image  proce.ssing,  to  multiply  first  the  pixel  values  such  that  the 
maximum  value  becomes  4095.  The  best  results  are  acquired  by  optimizing  the  exposure  during 
the  measuring  session  such  that  the  maximum  pixel  value  is  4095  in  the  original  picture.  Median 
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filtering  can  be  applied  in  order  to  correct  defective  or  hot  pixels.  Without  the  mentioned 
preprocessing  methods,  the  FEL-algorithm  still  resulted  in  a  significant  improvement.  Therefore, 
the  preprocessing  is  not  carried  out  in  most  of  the  cases. 

After  processing,  the  images  need  to  be  transformed  according  a  relation  given  by  <2.1>,  to  add  a 
gamma  correction  for  a  proper  monitor  display.  This  is  necessary  because  our  CCD-camera  is  a 
linear  recorder  of  the  scene  luminance.  In  most  ca.ses,  a  gamma  of  0.65  was  applied  with  good 
results.  When  a  larger  gamma  is  used,  then  the  bright  part  of  the  image  is  enhanced  at  the  cost  of 
the  dim  part  and  when  a  smaller  gamma  is  used,  the  other  way  around.  For  a  reliable  comparison 
of  the  original  image  with  the  processed  image,  the  original  image  had  to  be  displayed  with  a 
gamma  of  0.4,  which  is  a  proper  gamma  correction  for  most  of  the  monitors.  The  processing  with 
the  FEL-algorithm  was  carried  out  with  a  3x3  moving  average  filter,  w'hich  gave  the  best  results 
for  all  the  processed  images.  The  parameters  yj  and  k(,u,  (formula  <2.1.>)  have  been  calculated  so. 
that  the  dynamic  range  in  the  output  image  (without  contrast-enhancement)  runs  from  8  to  224.  A 
sufficient  large  contrast-enhancement  showed  to  be  possible  then  in  all  grey  level  ranges  without 
surpassing  the  range  from  0  to  255  in  the  final  output  image.  The  contra;  ts  have  been  enhanced 
till  no  further  improvements  could  be  observed  and  noi.se  or  artefact  phenomena  started  to 
increase.  A  value  for  '4=25  .showed  to  be  a  good  choice  for  most  of  the  12-bit  images.  Exceptions 
are  the  scenes  with  mainly  natural  vegetation  and  recorded  under  bright  sky  conditions  without 
photopic  filter. 

For  those  scenes  a  7^=15  was  a  good  choice.  With  the  parameter  v^,  the  contrast  enhancement  is 
adapted  such  that  no  artefacts  will  arise.  It  has  been  mentioned  already  in  chapter  two,  that  the 
parameter  v^  defines  a  threshold  in  the  local  variance  range.  Beyond  this  threshold  the  contrast 
enhancement  is  gradually  decreasing  to  one,  according  a  relationship  with  v„  as  a  variable. 

It  means  that  with  a  different  v^,  also  the  adaptation  relation  will  be  different.  Therefore,  some 
margin  in  the  values  of  y^  and  v„  will  exist,  also  in  case  of  optimal  processing.  Artefacts  first 
arise  in  the  pictures  at  the  horizon  and  they  often  can  be  observed  in  the  histogram  of  the 
processed  image  as  a  peak  at  greylevel  value  255.  When  the  contrasts  are  strongly  enhanced, 
other  artefacts  become  visible,  which  are  supposed  to  be  related  to  the  array  structure  of  the  CCD- 
sensor.  It  could  be  of  interest  to  investigate  these  artefacts  as  well  as  the  sources  of  noise  in  more 
detail  to  determine  the  limits  of  a  maximum  contrast  enhancement.  It  is  stressed  here  once  more 
that  the  contra.sts  in  the  images  are  affected  by  the  dynamic  range  reduction  due  to  the 
transformation  from  1 2  to  8  bit  (for  presentation  purposes  on  a  monitor). 
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The  results  obtained  with  the  FEL-algorithm  will  be  discussed  now. 

This  FEL-algorithm  gave  good  results  for  all  the  processed  images.  Images  with  a  dynamic  range 
more  that  8  bits,  recordings  with  facing  light  and  scenes  with  much  vegetation  were  hardly 
available  and  should  be  recorded  in  the  future.  It  was  already  mentioned,  that  a  special  effect  of 
the  shutter  might  arise  in  scenes  observed  through  haze.  Also  the  dynamic  range  transformation 
showed  to  be  not  optimal  for  images  with  a  very  limited  grey  level  range.  Due  to  the  structure  of 
the  algorithm  it  is  difficult  to  recognize  the  effects  of  the  various  parameters.  The  local  variance  is 
used  to  determine  the  adaptive  contrast  enhancement;  it  would  be  easier  to  use  the  local  difference 
in  the  logarithmic  domain,  which  is  equal  to  the  contrast  in  the  input  image.  The  way  in  which  the 
contrast  enhancement  is  realized  for  large  contrasts  depends  on  the  value  of  v^  and  the  local 
variance  in  the  logarithmic  domain.  Also,  this  makes  it  difficult  to  estimate  the  effect  of  v^.  but  at 
the  same  time  it  allows  some  margin  in  the  values  of  v^,  and  Yc-  The  FEL-algorithm  has  been 
compared  with  wellknown  image  processing  tools  like  histogram  equalizing  and  histogram 
stretching,  For  an  image  with  a  large  dynamic  range  (image  nr  db891005020l.  .see  Appendix  A), 
these  processing  techniques  have  been  applied  to  the  full  image  and  to  a  part  of  the  image;  in  the 
selected  part  only  a  limited  grey  level  range  was  pre.sent.  The  results  are  shown  in  respectively  the 
figures  5.1  \  for  the  full  image  and  in  the  figures  5.2x  for  the  lower  right  quarter  of  the  same 
picture.  Figure  5.1a  gives  the  original  image,  without  processing;  figure  5.1b  gives  the  result  of 
the  histogram  equalization  and  figure  5.1c  gives  the  result  with  the  FEL  image  processing 
algorithm. 

In  figure  5.2a  the  result  with  histogram  .stretching  is  given  for  the  lower  right  quarter  of  the  same 
picture  as  given  in  figure  5.1.  In  figure  5.2b  the  result  is  given  with  the  FEL  image  processing 
algorithm  for  the  same  part. 


Fij:ure  5.1a;  Original  12-bit  picture  without  processing. 


Figure  5.1b:  Histogram  equalization  applied  to  the  picture,  given  in  figure  5.1a. 


Figure  5.1c:  FEL  image  processing  algorithm  applied  to  the  picture,  given  in  figure  5.  la. 


Figure  5.2a:  Histogram  stretching  applied  to  the  lower  right  quarter  of  the  picture,  given  in  figure  5. 1 
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Figure  5.2b:  FEL  image  processing  algorithm  applied  to  the  lower  right  quarter  of  the  picture,  given  in 
figure  5.1a. 

With  histogram  equalization,  a  part  of  the  image  is  enhanced  at  the  cost  of  another  part  (see  fig. 
.5.1a  and  fig.  5.1b).  Details  in  parts  with  less  frequently  occurring  pixel  values  can  get  lost,  due  to 
histogram  equalization. 

The  contrast  stretch  function  can  only  be  used  for  images  with  a  limited  range  of  succeeding  grey 
values.  This  also  means  that  there  must  be  no  defective  photosites  causing  black  or  hot  pixels  or 
small  pixel  clusters.  The  contrast  stretch  function  is  not  appropriate  for  detail  contrast 
enhancement.  Even  when  applied  to  a  part  of  the  image  only  a  small  improvement  can  be 
ob.served  and  also  now  the  FEL  image  processing  algorithm  results  in  a  much  better  detailed 
picture  (see  fig.  5.2).  Comparing  figures  5.1c  and  5.2b  shows  that  the  FEL-algorithm  gives 
already  the  same  result  in  processing  the  complete  picture  as  when  it  is  applied  to  a  part  of  the 
image,  although  a  larger  range  is  available  then  for  the  contrast  enhancement  (this  is  to  be  seen  on 
the  monitor,  not  really  good  on  the  photographs).  It  can  be  concluded  that  an  algorithm  which 
makes  u.se  of  local  statistics  is  more  suited  to  enhance  an  entire  image  then  simple  histogram 
modification  techniques. 
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Because  the  pictures  in  figure  5.2.  are  zoomed  in,  some  of  the  sensor  array  structure  get  to  be 
seen. 

Now  results  with  examples  will  be  given  for  the  PEL  image  processing,  applied  to  images 
recorded  under  various  atmospheric  conditions. 

Figure  5.3  shows  an  image  that  was  recorded  with  a  photopic  filter  under  conditions  of  heavy 
atmospheric  haze.  The  visibility  varied  in  the  vertical  direction  from  1.3  to  about  4.5  km.  Figure 
5.4  shows  an  image  recorded  at  about  the  same  time  and  with  about  the  same  visibility  (1.4  to 
about  4  km),  but  now  recorded  without  the  photopic  filter.  It  can  be  seen  that  after  processing  a 
significantly  better  result  is  obtained,  especially  in  regions  in  which  larger  contrasts  are  expected 
with  near  infrared  radiation  (horizon  and  regions  of  vegetation).  Without  the  photopic  filter  the 
near  infrared  radiation  contributes  to  better  contrasts  in  regions  of  vegetation  because  of  a  larger 
reflection  in  this  spectral  region.  Also  the  better  transmission  of  the  near  infrared  wavelengths  in 
dense  fog  is  of  importance,  especially  for  the  larger  distances.  In  general,  a  better  result  is 
obtained  without  a  photopic  filter  for  visibilities  below  8  km,  except  in  sunlit  scenes.  After 
processing,  small  rectangular  areas  become  visible  due  to  the  large  enhancement.  These  kinds  of 
artefacts  and  noise  effects  are  more  clearly  visible  on  the  monitor.  Figure  5.5  shows  an  image, 
taken  with  a  photopic  filter  at  a  visibility  of  7  km.  It  can  be  seen  that  this  image  should  have  been 
recorded  without  photopic  filter.  Figure  5.6  shows  an  image  recorded  with  slight  facing  light'  at  a 
visibility  of  30  km.  It  can  be  seen  that  the  12-bit  camera  system,  provided  with  image  processing, 
.still  can  make  perceivable  more  details  in  this  case.  Moreover,  it  results  in  a  sharper  picture  with  a 
more  comfortable  vision.  Figure  5.7  shows  a  picture,  taken  during  rain  at  a  visibility  of  about  10 
km.  After  processing  some  optical  distortions  due  to  the  rain  become  perceivable.  Figure  5.8 
shows  a  picture  with  mainly  vegetation,  taken  with  photopic  filter  at  a  visibility  of  about  17  km.  It 
is  an  example  of  processing  with  a  lower  value  for  the  contra.st  enhancement  parameter  (Yc='  -5)  in 
order  to  prevent  artefacts.  In  .scenes  with  much  vegetation,  artefacts  sooner  arise  then  in  the  city 
views'.  The  picture  in  figure  5.8  gives  also  a  good  example  of  processing,  which  results  in  an 
extended  visual  range  and  a  remarkable  increased  comfortable  vision.  Figure  5.9  finally  shows 
what  happens  when  the  photopic  filter  is  not  used  in  a  sunlit  scene  under  conditions  of  good 
visibility. 
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5.4a 


5.4b 


Figure  5.4: 


Picture  taken  without  photopic  filter;  visibility  from  1.4  to  4km  (image  dr88 1 1070403. 
appendix  A).  Top  before  and  bottom  after  processing. 
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5.5b 


Figure  5.5:  Picture  taken  with  photopic  filter;  visibility  7  km  (image  m8901 270201,  appendix  A).  Top 
before  and  bottom  after  processing. 
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5.6a 


Picture  taken  with  photopic  filter  at  a  very  good  visibility  of  30  km  (image  cc88101(X)60l 
appendix  A).  Top  before  and  bottom  after  processing. 


Figure  5.6: 
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5.6b 
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5.7a 


5.7b 


Figure  5.7: 


Picture  taken  during  rain;  visibility  10  km  (image  st8809230101,  appendix  A).  Top  before 
and  bottom  after  processing. 
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5.8a 


Figure  5.8: 


Picture  with  mainly  vegetation,  taken  with  photopic  filter;  visibility  17  km  (image 
dr881 1070501,  appendix  A).  Top  before  and  bottom  after  processing. 
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Figure  5.9;  Unprocessed  picture,  taken  without  photopic  filter  at  bright  weather;  visibility  30  km  (image 
sp8806060l03,  appendix  A). 

5.2  Analysis  of  the  image  file 

Some  relations  between  parameters  of  the  FEL-algorithm  and  meteorological  parameters  have 
been  studied,  in  view  of  possible  automation  of  the  image  processing  parameters.  It  has  been 
mentioned  before  that  the  parameters  and  kgu,  can  be  calculated  from  the  minimum  of  the 
histogram,  y^.  can  be  kept  constant  (at  first  instance)  and  the  adaptation  parameter  v„  can  be 
determined  from  the  calculated  local  variance.  This  local  variance  depends  on  the  atmospheric 
conditions  and  the  scene.  The  minimum  in  the  histogram  distribution  of  the  grey  values  depends 
to  a  great  extent  on  the  scattering  of  light  in  the  atmosphere  (see  for  the  contribution  of  the 
background  luminance  formula  4.2).  In  general,  this  minimum  value  will  increase  with  increasing 
Relative  Humidity  and  with  decreasing  Visibility.  From  studies  in  the  field  of  atmospheric 
scattering  it  is  known  that  the  parameter  Relative  Humidity  /  Visibility  correlates  well  with 
atmospheric  scattering,  especially  forward  scattering  [8]. 

Some  image  parameters  will  be  correlated  now  with  the  scattering  parameter  Wsc,  which 
therefore  is  defined  as  Wsc  =  Relative  Humidity  /  Visibility. 


Min  greylevel  (12bit  AD  units) 
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Figure  5.10  shows  the  minima  of  the  grey  level  distributions  as  a  function  of  this  scattering 
parameter  Wsc  for  all  the  processed  images.  The  given  straight  line  is  obtained  by  a  first  order 
linear  regression  of  the  logarithmic  values  of  the  relevant  parameters.  A  correlation  coefficient  of 
about  S9%  was  found.  Table  5.1  gives  the  slopes  and  interceptions  of  such  regression  lines  with 
the  correlation  coefficients  for  various  selections. 


10000 
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1 

Figure  5.10:  Minimum  value  in  the  greylevel  distributions  as  function  of  the  scattering  parameter  Wsc. 
for  all  the  processed  pictures. 
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Table  5.1:  Relations  between  atmospheric  conditions  and  some  image  parameters. 


1 

2 

3 

4 

5 

6 

7 

8 

nr 

parameter 

range 

filter 

a 

b 

correlation 

number  pict. 

1 

minimum 

0-30 

0/1 

0.5 

2.3 

0.89 

173 

2 

0-30 

0 

0.4 

2.5 

0.89 

46 

3 

0-30 

1 

0.6 

2.2 

0.9 

127 

4 

0-7 

0/1 

0.4 

2.5 

0.85 

107 

5 

0-7 

0 

0.3 

2.6 

0.86 

79 

6 

8-30 

0/1 

0.7 

2.1 

0.7 

66 

7 

loc.var 

8-30 

1 

0.7 

2.1 

0.73 

50 

8 

vn/957f 

0-30 

0/1 

-1.1 

-2.0 

-0.79 

143 

9 

vn/957f 

0-30 

0 

-0.9 

-2.0 

-0.87 

37 

10 

vn/95% 

0-30 

1 

-1.2 

-1.9 

-0.79 

no 

n 

vn/chsn. 

0-30 

0/1 

-0.6 

-3.0 

-0.57 

159 

12 

vn/chsn 

0-30 

0/1 

-0.7 

-3.0 

-0.71 

88 

13 

Vn/chsn 

0-30 

1 

-0.7 

-3.0 

-0.72 

68 

14 

Vn/chsn 

0-30 

1 

-0.6 

-3.0 

-0.65 

36 

15 

Vn/chsn 

0-7 

0/1 

-0.6 

-3.1 

-0.54 

52 

16 

Vn/chsn 

8-30 

0/1 

-0.8 

-2.8 

-0.59 

36 

In  table  5.1  the  respective  columns  represent: 


1 .  nr:  is  the  graph  number  for  reference  purpo.ses. 

2.  parameter:  is  the  histogram  parameter,  which  is  correlated  with  the  scatter  parameter  Wsc. 
The  local  variance  parameter  v„/95%  is  the  value  in  the  local  variance  histogram  distribution 
at  the  957r  point.  The  local  variance  parameter  v„/chsn  is  the  value  of  the  local  variance, 
which  is  actually  used  in  the  image  processing. 

.3  range:  is  the  visibility  range  (in  km),  which  is  taken  into  account  for  the  given  picture 
selection. 

4  filter:  gives  the  picture  selection  with  respect  to  the  used  filter;  0  are  pictures  taken  without 
filter  and  1  pictures  taken  with  a  photopic  filter. 

5  a:  gives  the  slope  of  the  first  order  regression  line  with  the  parameters  on  a  logarithmic  scale. 

6  b:  gives  the  intersection  at  the  vertical  axis  (with  the  histogram  parameter)  of  the  first  order 
regression  line,  the  parameters  are  supposed  to  be  on  a  logarithmic  scale. 

7  correlation:  correlation  coefficient  obtained  with  first  order  regression. 

8  number  pict.:  number  of  pictures  selected  for  the  respective  graphs. 
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The  graph  numbers  i  to  7  give  the  regression  line  parameters  for  the  graphs  with  the  minimum 
grey  levels  in  the  histogram  distributions  (graph  nr  1  is  pictured  in  figure  5.10);  the  nrs  8  to  10 
give  these  parameters  for  the  local  variance  parameter  v„  at  95%  of  its  histogram  distributions; 
the  nrs  11  to  16  for  the  local  variance  parameter  v„,  which  has  been  actually  used  in  the 
processing.  The  actually  used  values  of  v^  at  the  corresponding  percentages  in  their  histogram 
distributions  can  be  found  in  column  1 1  of  the  list  of  processed  pictures  in  appendix  A. 

In  graph  nr  12  of  table  5.1  only  the  pictures  with  a  contrast  multiplier  Yc=25  are  considered.  In 
graph  nr  1 3  only  the  pictures  concerning  the  city  views  with  a  contrast  multiplier  Yc=25  have  been 
taken  into  account.  In  graph  nr  14  the  same  pictures  as  in  nr  12  have  been  considered,  but  pictures 
taken  with  the  wind  directions  west  and  southwest  have  been  omitted.  In  the  nrs  15  and  16  only 
the  pictures  with  a  contrast  multiplier  Yc=25  have  been  taken  into  account. 

The  visibilities  (especially  of  the  'city  view.s’)  have  been  determined  by  eye,  in  most  cases.  This 
visibility  was  supposed  to  be  equal  to  the  distance  to  a  relevant  reference  point,  which  just  could 
be  seen  by  a  trained  observer.  For  visibilities  equal  to  or  better  than  30  km  no  such  references 
were  available. 

The  nearest  objects  in  the  pictured  scenes  are  at  a  distance  of  about  1  km  from  the  camera,  which 
always  results  in  an  offset  of  the  grey  level  range.  This  offset  results  from  a  contribution  of  light 
caused  by  atmospheric  scattering  (see  also  formula  4.2).  The  lowest  minimum  grey  level  that  has 
been  recorded,  amounted  to  a  value  80  (out  of  4095  in  12-bit  AD  units)  and  the  highest  about 
2800  (see  figure  5.10).  It  might  be  clear  that  the  dynamic  range,  in  general,  depends  on  the  scene 
itself.  In  our  recordings  the  maximum  grey  level  value  always  originated  from  the  sky  beyond  the 
horizon  and  the  minimum  value  from  the  vegetation  on  the  foreground.  The  dynamic  range 
parameters  Yd  and  k^^  depend  on  the  ratio  between  the  maximum  and  the  minimum  luminance,  so 
in  our  case  no  explicit  dependence  on  the  scene  may  be  expected. 

From  table  5.1  it  can  be  seen  that  the  correlation  with  the  scatter  parameter  Wsc  is  better  for  the 
minimum  greylevel  than  for  the  adaptation  parameter  Vp 

In  figure  5.11  the  local  variance,  at  the  95%  points  of  the  histogram  distributions,  is  given  as 
function  of  the  scatter  parameter  Wsc  (Rel.HumA'is)  for  all  the  processed  pictures  (graph  nr  8  in 
table  5.1). 


local  variance  (95%  histogram) 
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local  variance  adaptation  threshold 


Figure  5.11:  Local  variance  at  95%  point  of  the  distributions  as  function  of  the  scatter  parameter  Wsc.  for 
all  the  processed  pictures. 

The  selections  of  pictures  taken  with  photopic  filter  have,  in  general,  a  lower  value  for  the 
intersection  parameter  'b'  and  a  steeper  slope  'a'  than  the  pictures  taken  without  filter. 

The  correlation  of  the  adaptation  parameter  v„  with  the  scatter  parameter  Wsc  is  substantially 
better  for  a  value  at  a  fixed  percentage  (95%)  of  its  histogram  distribution  than  for  the  actually 
used  values,  which  gave  (in  a  subjective  sense)  the  best  processing  results.  The  local  variance  is 
related  to  the  local  contrast,  which  is  determined  by  the  contrasts  of  the  objects  in  the  scene  and 
by  the  contrast  loss  in  the  atmosphere.  The  scattering  increases  when  the  Relative  Humidity 
increases  and  the  Visibility  decreases.  The  parameter  Wsc  =  Relative  Humidity  /  Visibility 
therefore  is  a  relevant  choice  to  correlate  with.  It  can  be  seen  from  graph  1 1  in  table  5.1  that, 
when  the  actually  used  v^  is  related  to  this  scatter  parameter  Wsc  for  all  the  images,  there  is 
hardly  a  relevant  correlation  (R=0.57).  It  must  be  emphasized  here  that,  when  an  image  is 
processed  with  a  different  or  different  dynamic  range  reduction  (output  range),  the  value  of  v^, 
will  also  be  different.  If  only  the  images  with  a  Y(.“25  are  taken  into  account  and  the  same 
dynamic  range  reduction  is  applied,  then  the  correlation  is  better  (see  figure  5.1 1.  graph  nr  12, 
R=0.7I). 
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Because  the  images  are  recorded  near  the  coast  (North  Sea),  the  scattering  aerosol  particles  are  of 
mainly  maritime  nature,  especially  when  the  wind  direction  is  between  south  and  west.  Omitting 
all  the  images  recorded  during  such  a  wind  direction,  results  in  a  lower  correlation  coefficient 
(R=0.65),  however.  This  result  may  be  due  to  the  very  few  points  left  in  this  case.  Furthermore, 
there  is  a  fluctuation  to  expect  in  the  parameter  'chosen  v^',  because  the  perception  of  the  best 
processing  quality  is  subjective  by  nature.  The  local  variance  at  95  %  of  its  histogram 
distributions  (see  chapter  two)  gives  a  better  correlation. 

It  might  be  concluded  that  the  correlation  of  the  minima  in  the  'grey  level'  histogram  distributions 
(and  to  a  lesser  extent  the  parameter  v^)  with  the  scatter  parameter  Wsc  is  sufficient  to  use  it  for  a 
parameter  setting  of  the  image  processing.  However,  it  seems  more  reliable  to  deduce  such 
parameters  directly  from  the  histogram  distribution  of  the  grey  levels,  if  such  histograms  are 
available. 

5.3  Visibility  range  improvement 

It  was  already  mentioned  in  chapter  four  that  it  is  very  hard  to  determine  the  visibility  range 
improvement  due  to  image  processing.  A  simple  comparison  of  a  monitor  presentation  of  natural 
scenes  before  and  after  the  image  processing  poses  special  requirements  on  the  composition  of  the 
scene.  Only  a  few  12-bit  images  could  be  used  for  this  purpose.  The  distance  to  only  few  objects 
in  the  scene  are  known.  One  must  find  an  object  at  known  distance  in  the  processed  image  seen 
against  the  horizon  and  a  similar  object  in  the  original  image  with  about  the  same  contrast.  The 
apparent  visibility  range  improvement  due  to  image  processing  can  be  estimated  then.  For  the 
contrast  improvement  due  to  processing,  it  is  necessary  to  choose  an  object  in  the  processed 
image,  because  an  object  that  becomes  perceivable  only  after  processing  has  the  maximum 
contrast  improvement,  while  the  other  objects  can  be  processed  in  an  adaptive  way.  The  results 
for  a  few  images  are  given  in  table  5.2.  The  visibility  used  in  table  5.2  is  the  visibility  estimated 
by  eye  during  the  recordings;  R1  and  R2  are  the  distances  to  the  objects  in  respectively  the 
images  before  processing  and  after  processing;  Cr  and  Cr'  the  contrasts  of  the  chosen  objects  in 
respectively  the  images  before  processing  and  after  processing.  These  contrasts  are  determined  by 
averaging  the  grey  levels  in  and  around  the  chosen  objects.  The  'visibility  range  improvement'  VI 
can  be  calculated  then  with  Cr'/Cr  from  formula  4.6,  assuming  that  Co=-l  (a  complete  black 
object  according  formula  4.1).  The  dependence  of  the  'visibility  range  improvement'  on  the  initial 
contrast  Co  is  not  analysed. 
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It  can  be  seen  that  the  'visibility  range  improvement’  due  to  the  image  processing  varies  from  1 .2 
to  3.5  for  the  given  examples  and  premises.  The  observed  visibility  (by  eye)  of  course  lies  in 
between  the  distance  to  the  chosen  object  before  processing  (Rl)  and  the  one  after  processing 
(R2). 


Table  5.2;  Visibility  range  improvement  due  to  image  processing. 


filename 

visi¬ 

bility 

Rl 

R2 

R2/R1 

Cr 

Cr 

Cr'/Cr 

VI 

St  8809210401 

3.4 

>1.2 

3.4 

<2.8 

-0.03 

-0.08 

2.7 

1.4 

dr  8811090101 

4.0 

3.4 

>4.2 

>1.2 

-0.02 

-0.13 

7.0 

2.0 

dr  8811090102 

4.0 

3.4 

>4.2 

>1.2 

db 89072601 01 

2.5 

>1.2 

<3.4 

<2.8 

-0.10 

-0.25 

2.5 

1.7 

db  8907260201 

4.5 

3.4 

4.2 

>1.2 

-0.03 

-0.19 

6.3 

2.1 

cc  8810110501 

2.7 

3.4 

>4.2 

>1.2 

-0.01 

-0.04 

4.0 

1.4 

cc  8810110503 

3.4 

>3.4 

>4.2 

>1.2 

-0.03 

-0.16 

5.3 

1.9 

cc  8810110504 

3.5 

>3.4 

>4.2 

>1.2 

-0.03 

-0.17 

5.7 

2.0 

cc  8810110602 

4.2 

3.4 

>4.2 

>1.2 

-0.01 

-0.10 

10.0 

2.0 

dr  8811070601 

4.0 

1.2 

4.2 

3.5 

-0.01 

-0.09 

9.0 

1.9 

dr  8811070501 

1.3 

1.2 

4.2 

3.5 

-0.02 

-0.07 

3.5 

1.5 

dr  8811070703 

3.6 

3.4 

>4.2 

>1.2 

-0.02 

-0.07 

3.5 

1.5 

dr  8811070403 

1.4 

1.2 

4.2 

3.5 

-0.03 

-0.17 

5.7 

2.0 

sp  8810040401 

4.0 

>1.2 

3.5 

<2.9 

-0.01 

-0.04 

4.0 

1.4 

St  8809220201 

2.8 

0.5 

1.2 

2.4 

-0.03 

-0.08 

2.7 

1.4 

St  8809220501 

4.5 

1.2 

4.2 

3.5 

-0.03 

-0.15 

5.0 

1.8 

St  8809220401 

3.0 

<3.5 

>3.5 

>1.0 

-0.1 

-0.18 

1.8 

1.3 

St  8809210101 

3.4 

1.2 

3.4 

2.8 

-0.04 

-0.09 

2.3 

1.3 

The  accuracy  of  the  given  values  is  rather  low,  because  often  different  objects  have  been 
compared  before  and  after  processing.  Nevertheless,  they  give  an  indication  of  what  improvement 
may  be  expected.  In  chapter  four  a  better  experiment  with  black  and  white  screens  is  already 
proposed.  Recordings  could  already  be  made  from  black  and  white  screens  with  increasing 
dimensions  at  increasing  distances,  for  one  type  of  visibility.  These  screens  with  an  initial  contrast 
Co=0.89  were  located  at  distances  of  130,  170  and  290  metres  from  the  camera.  The  dimensions 
were  such,  that  in  the  pictures,  the  screens  all  matched  with  an  equal  number  of  pixels  on  the 
CCD-sensor.  The  results  are  given  in  table  5.3.  TTie  'visibility  range  improvement'  VI  (according 
formula  4.6),  is  in  the  same  range  as  the  former  estimated  improvement  from  the  monitor 
presentation,  notwithstanding  the  screens  were  too  small  to  get  an  accurate  contrast  value  by 
averaging  the  grey  levels. 
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Table  5.3:  Visibility  range  improvement  from  black  and  white  screens. 


filename 

screen 

Cr 

Cr' 

VI 

R2/R1 

db8904130101 

2 

-0.06 

-0.3 

2.3 

>1.7 

db8904 130203 

2 

-0.06 

-0.3 

2.3 

2.2 

db8904 130301 

3 

-0.06 

-0.3 

2.3 

>1.3 

The  results  from  table  5  2.  and  5.3.  show  that  the  'visibility  range  improvement'  lies  in  between  1 
and  3.  Table  5.3  represents  only  one  weather  type:  fog  with  a  decreasing  density. 

5.4  Processing  procedure 

A  straightforward  procedure  for  obtaining  optimum  processing  can  be  given  now  for  the 
algorithm,  evaluated  in  this  study. 

The  calculation  of  the  local  variance  is  necessary  for  the  adaptive  contrast  enhancement.  This 
makes  it  necessary  to  store  two  files,  one  file  containing  the  local  mean  and  one  file  containing 
the  local  variance  of  each  pixel.  The  parameter  showed  to  be  rather  constant  for  the  here 
considered  kind  of  scenes  (see  appendix  A:  list  of  processed  images  and  processing  parameters). 
The  start  values  of  some  parameters  may  be  deduced  from  the  relative  humidity/visibility  (see  the 
correlation  and  regression  line  parameters  in  table  5. 1 ). 

A  straight  forward  procedure  now  is: 

The  range  transform  parameters  and  k^y,  can  be  determined  from  the  input  range  (min- 
max  range  of  the  grey  level  hi.stogram  distribution)  and  the  output  range  (8-224). 

Calculate  the  local  difference  file  and  local  variance  file  with  a  very  low  value  of  v^  (eg 
lE-7). 

Take  as  start  value  for  the  adaptation  parameter  v^  the  value  at  the  959c  point  in  the  local 
variance  histogram  distribution. 

Lower  the  value  of  till  all  artefacts  disappear,  while  applying  a  3x3  window  and  a 
maximum  value  for  the  contrast  multiplier  (y(.=25). 

Eventually,  optimize  the  value  for  the  contrast  multiplier. 
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6  CONCLUSIONS 

The  optimum  contrast  enhancement  factor  is  rather  constant  for  the  considered  scenes. 
Extending  the  database  with  pictures  from  other  scenes,  especially  with  large  dynamic  range, 
is  recommended. 

For  pictures  with  a  very  narrow  histogram  of  the  greylevel  distribution,  the  range 
transformation  can  better  be  performed  in  advance  of  the  local  contrast  enhancement. 

In  case  that  details  in  vegetation  parts  of  pictures  are  important,  the  use  of  a  photopic  filter  is 
necessary  at  visibilities  beyond  about  7  km. 

The  perception  of  visibility  range  improvement'  by  image  processing  depends  among  others 
on  the  scene,  the  weather  type  and  the  observer.  This  extension  of  the  visibility  range  due  to 
processing  is  expected  to  lie  between  1 .5  and  3  times,  at  least  for  weather  situations  with  a 
considerable  loss  of  contrast  in  the  atmosphere. 

It  is  recommended  to  quantify  the  visibility  range  extension.  This  may  be  possible  in  a 
special  experiment  with  a  large  number  of  well  defined  targets  at  different  distances.  It  is 
important  to  quantify  the  usefulness  of  image  processing  under  various  atmospheric  and 
battle  field  conditions. 

It  is  recommended  to  study  the  positive  effects  of  contrast  enhancement  on  the  probability  of 
detection  and  recognition  of  military  objects. 


TNO  report 


Page 

45 


7  EPILOGUE 

Many  of  the  results  and  insights  of  the  study,  described  in  this  report,  have  been  used  for 
demonstrating  the  power  of  the  image  processing  algorithm  in  a  number  of  other  CCD-imaging 
projects. 

In  order  to  be  able  to  add  pictures  from  various  scenes  to  the  image  database,  a  PC-controlled 
camera  system  has  been  built  with  a  modular,  thermo-electric  cooled  camera  head  [9,10,1 1].  A 
large  number  of  pictures  have  been  taken  with  this  camera  system,  among  others,  during  the 
battlefield  trial  BEST  TWO  at  Mourmelon,  France  [10]. 

'J 

All  the  pictures  dealt  with  in  this  report  have  been  processed  on  a  VAX  computer  and  a  1“S  image 
proce.ssing  station.  For  most  of  the  interesting  applications,  a  stand  alone  system  with  real  time 
processing  will  be  demanded.  The  processing  time  with  the  described  algorithm  [2],  is  a  serious 
drawback  in  eg  a  PC  controlled  stand  alone  station.  A  general  purpose  parallel  computing  system 
was  used  for  the  development  of  a  real  time  processing  system  [15].  With  four  transputers  a 
processing  time  of  about  1  second  per  frame  was  realized.  Most  of  the  Mourmelon'  images  have 
been  processed  with  this  parallel  system  [10].  Variations  of  the  algorithm  have  been  studied  with 
real  time  hardware  implementation  in  mind  [12,13].  With  fast  electronic  components  a  real-time 
video  processor  has  been  realized,  which  accepts  any  standard  analog  video  signal  [14].  A  global 
range  transformation  is  implemented  in  this  unit  directly  on  the  analog  video  signal  with 
programmable  gain  and  offset.  After  digitization  local  adaptive  contrast  enhancement  is  carried 
out. 

Finally,  within  the  framework  of  a  national  technology  project,  a  12-bit  real  time  processing 
system  has  been  developed  and  realized,  together  with  a  number  of  real  time  filtering  modules 
[15].  This  complete  processing  system  (with  filtering  and  averaging  options)  is  part  of  a 
sophisticated  EBCCD  camera  system,  but  it  can  accept  also  any  standard  analog  video  signal. 
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Legend  of  list  of  processed  recordings. 

1-  filename:  abyymmddnncc  (eg  dr881 1030101) 
with: 

ab=tape  /  yy=year  /  mm=month  /  dd=day  /  nn=picture  nr/  cc=code 

code:  01=with  photopic  filter,  02=id  with  8  pictures  summed,  03=without  filter  /  04=id  with 
8  pictures  summed,  05=dark  picture. 

2-  lum:  Illuminance  in  kilolux 

3-  zi:  Visibility  in  km 

4-  hu:  Relative  Humidity  in  % 

5-  T :  Temperature  in  degrees  Celsius 

6-  wo:  cloudiness 

7-  wt:  weather  type 

with: 

1  =  rain 

2  =  clear  with  some  clouds 

3  =  sunlit  scene 

4  =  drizzle 

5  =  fog 

6  =  fog,  decreasing  density 

7  =  haze  at  some  distance 

8  =  some  local  cloudiness 

9  =  haze 

10  =  decreasing  haze 

1 1  =  clear 

1 2  =  dark  weather 

8-  Yj  exponent  of  range  transformation 

9-  k:  constant  in  range  transformation 

10-  contrast  multiplier  in  logarithmic  domain 

1 1-  vn/%:  local  variance  in  logarithmic  domain  for  the  given  Vr  in  its  histogram  distribution; 
this  value  has  actually  been  used  in  the  processing. 

12-  min:  lowest  grey  level  value 

13-  vn/95:  local  variance  in  logarithmic  domain  for  the  95‘7c  point  in  its  histogram  distribution. 

14-  S:  description  of  the  scenes 
with: 

1  city  view  1  with  Bronovo  hospital 

2  cityview  2  with  barrack  building 

3  cityview  3  with  central  station 

4  cityview  4  with  Esso  building 

5  white  house  between  dunes/woods 

6  small  building  in  dunes  near  TNO-FEL  parking 

7  read  roof  /  a-hove 

8  parking  STC 

9  roof  STC 
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15- 


opm:  remarks 
with: 

1  =  visibility  measured  by  AEG  point  visibility  meter  near  TNO-FEL  laboratory 

2  =  not  used  in  the  statistics 

3  =  bad  part  in  picture 

4  =  white  vegetation  due  to  omission  of  photopic  filter 

5  =  application  of  NIR  passing  filter  RG695 

6  =  no  information  perceivable  after  processing 

7  =  not  standard  output  window  (which  is  8,224) 

8  =  grey  level  scaling  before  processing 

9  =  no  homogeneous  visibility  in  vertical  picture  direction 

10  =  facing  light  present  in  picture 

1 1  =  7*7  filter  window  applied  in  local  contrast  enhancement  (standard  is  3*3) 

1 2  =  complete  picture  below  the  horizon 
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f ilenaam 

M 

hu 

1 

■ 

I 

i 

7d 

k 

7c 

vn/  Z 

min 

vn/95 

s 

opm 

rn8812020304 

1.4 

0.2 

99 

-  2 

a 

5 

7  .  53 

*.13e-25 

40 

26-5/95 

2280 

2.26-5 

2 

rn6812020402 

1.4 

0.2 

99 

*2 

8 

5 

7 . 80 

1.63a-26 

7 

7 

2650 

2.7e-5 

2 

rn8812020603 

1  2 

0.2 

99 

-2 

8 

5 

8. 17 

2.87e-27 

25 

26-5/95 

2296 

2.0e-5 

1 

rn8812020501 

1 . 2 

0.2 

99 

-  2 

8 

5 

a .  02 

2.81e-27 

25 

36-5/95 

2650 

3  ,  06-5 

1 

- 

ec8810110101 

5.3 

0.2 

99 

5 

8 

5 

7 . 09 

5.33e-2A 

25 

3a-5/95 

2560 

2.8e-5 

2 

6 

CC8810110103 

3  1 

0.2 

99 

8 

1  .  00 

5.476-20 

25 

0.1 

2489 

- 

2/6 

st8e09220101 

0 . 2 

99 

10 

8 

5 

7 . 09 

5.33e-24 

7 

7 

2560 

2.9e-5 

1 

6 

rn8812020101 

2.1 

0 . 2 

99 

-2 

6 

5 

6 . 80 

6. 406-23 

35 

46-5/95 

2460 

4  .  Oe- 5 

2 

- 

rn8812020203 

1  .  4 

0.2 

99 

-  2 

8 

7 . 60 

2.266-25 

50 

3e-5/95 

2280 

2.66-5 

2 

- 

rn8901240103 

1 . 0 

0 . 5 

92 

0 

8 

1 . 24 

7.00e-3 

10 

le-3/85 

2172 

- 

8 

8 

rn89012^0201 

1 . 0 

0 . 5 

92 

0 

8 

1 . 24 

6 . 50e-3 

10 

96-4/87 

2449 

8 

S 

rn8902090303 

33 

0 . 5 

99 

0 

4 

5 

4  .  05 

5.21e-13 

35 

7e-5/85 

1799 

4 . 26-4 

9 

- 

rn890209Q401 

22 

0  .  5 

99 

0 

4 

5 

4  .  42 

2.41e-14 

30 

6e-5/85 

1927 

5  .  le-4 

9 

- 

db8904l30703 

5.5 

0.5 

98 

5 

8 

5 

1  .  00 

5.00e-2 

10 

56-2/- 

2120 

- 

2/8 

44 

0 . 7 

99 

0 

4 

5 

5  .  56 

1 . 906-18 

25 

46-5/95 

2248 

3.6e-5 

2 

- 

rn8902090603 

66 

0  .  7 

99 

1 

3 

5 

*  .  58 

8.62e-15 

25 

6e-5/95 

1847 

6  Oe-5 

2 

- 

sp88100i.0201 

22 

0.7 

98 

1  2 

7 

•> 

7 . 98 

3.386-27 

20 

26-5/85 

2697 

2 . 9e-5 

2 

- 

sp88100<.0202 

26 

0.  8 

98 

12 

7 

9 

1.726-24 

20 

26-5/90 

2535 

2 . 6e-5 

2 

- 

sp8810040101 

7 

0  .  8 

98 

12 

7 

9 

6.  52 

8.126-22 

25 

26-5/85 

2457 

3. Oe-5 

2 

- 

sp8810030101 

18 

0  .  8 

96 

1  0 

8 

6 

8  .  38 

1.16e-28 

20 

6e-5/99 

2752 

- 

2 

- 

rna90124030: 

0.9 

1 . 0 

92 

0 

8 

5 

1 . 27 

6. OOe-3 

10 

5e-4/85 

1940 

8 

8 

rn89012«0403 

0 . 9 

1  .  0 

92 

0 

8 

5 

1 . 2* 

7.00e-3 

10 

96-4/75 

1940 

8 

8 

CC8810110203 

7 

1  ,  2 

95 

5 

8 

10 

4  .  42 

4 .236-14 

25 

86-5/95 

1700 

8 . 3e-5 

2 

- 

sp8810G30201 

28 

1 . 2 

96 

1  0 

8 

6 

6.  73 

1 . 096-22 

20 

7e-5/99 

2496 

- 

2 

- 

sp881003030: 

28 

1 . 2 

96 

1  0 

8 

6 

6  .  15 

1.42e-20 

15 

2e-5/85 

2381 

- 

2 

- 

CC8810110201 

7 

1 .  3 

95 

5 

7 

10 

3  .  48 

6, 236-11 

25 

9e-5/95 

1570 

8  .  7e-5 

2 

- 

dr68101r.0103 

22 

1 .  3 

96 

9 

4 

6 

3  .  43 

2 . 40e-10 

25 

6e-5/5 

1375 

7 . 8e-4 

2 

1 

dr8811070501 

16 

1.3 

98 

2 

6 

3  .  40 

1 . 196-10 

25 

le-4/70 

1536 

- 

1 

9 

dr8811070403 

5,  5 

1  .  4 

98 

2 

2 

5 

3  .  18 

7 . 57e-10 

30 

le-3/60 

1434 

- 

1 

9 

dr8811070301 

7  .  5 

1  ,  4 

98 

2 

2 

5 

3  .  52 

4 . 49e-ll 

25 

le-5/70 

1587 

- 

2 

9 

dr8811070203 

5.5 

1.5 
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